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Abstract Modified gravity models have been constantly proposed with the pur¬ 
pose of evading some standard gravity shortcomings. Recently proposed by A.H. 
Chamseddine and V. Mukhanov, the Mimetic Gravity arises as an optimistic alter¬ 
native. Our purpose in this work is to derive Tolman-Oppenheimer-Volkoff equa¬ 
tions and solutions for such a gravity theory. We solve them numerically for quark 
star and neutron star cases. The results are carefully discussed. 

Keywords Modified gravity • compact stars ■ Tolman-Oppenheimer-Volkoff 
equations 


1 Introduction 

Since the late 90’s, it is known that our universe is undergoing a phase of accel¬ 
erated expansion. The first evidence came from Supernovae la observations m- 
Thenceforth, evidences from baryon acoustic oscillations [3] and cosmic microwave 
background anisotropies [3] seem to strengthen this cosmological scenario. If one 
assumes General Relativity (GR) as the background theory of gravity for a cos¬ 
mological model, there is a necessity of considering ~ 70% of the universe is filled 
with some exotic fluid named dark energy, which would be the responsible for 
the cosmic acceleration. However such a universe set up, although provide good 
agreement with observations, fails in explaining or describing some issues, like 
hierarchy, coincidence and cosmological constant problems [5ll6ll7[[8lCTlO] . among 
others. Some of those issues may, in principle, be solved (or at least evaded) by 
assuming alternative models of gravity. 

An important class of alternative models arise from the assumption of a gen¬ 
eralization of the Ricci scalar dependence on Einstein-Hilbert action, from which 
gravity field equations are derived, namely f(R) theories [llj-[13_i. 
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A generalization of f(R) theories has been recently proposed by T. Harko et al. 
m , for which the gravitational part of the action depends not only on a function 
of R but also on a function of T, the trace of the energy-momentum tensor. 

Those theories have been intensively tested in the cosmological level (check, for 
instance, m -[20]). However, alternative models of gravity should also be tested 
in the astrophysical level. In fact, strong gravitational fields found in relativistic 
stars could discriminate standard gravity from its generalizations. 

In what concerns f(R) theory astrophysical tests, in [3T], the authors have ob¬ 
tained f(R) = R+aR 2 Tolman-Oppenheimer-Volkoff (TOV) equations, which were 
solved numerically using adaptive Gaussian quadrature. A general f(R ) model 
with uniformly collapsing cloud of self-gravitating dust particles was studied in 
[22 :. Such an investigation could be used as a first approximation to the modifi¬ 
cation that stellar objects can suffer in modified theories of gravity. In [23], the 
gravitational field of stars were studied both analytically and numerically and 
generalized TOV equations have been derived. 

f(R,T) gravity has also been tested in such a level. In [23], the authors have 
analyzed the dynamics of a spherically symmetric collapsing star surrounding in 
locally anisotropic environment with expansion-free condition, while in [25] the 
authors have focused on the study of the stability of such objects. 

Recently, another kind of alternative gravity model has been proposed, for 
which the metric is not considered a fundamental quantity. Instead, it is taken 
as a function of an auxiliary metric g^ and a scalar field <j>, contemplating the 
so called Mimetic Gravity (MG) [26]. Note that such a dependence of the metric 
makes the application of the variational principle in the model action to yield 
more general equations of motion (EoM) than that of purely Einsteinian relativity 
theory. 

In |27| . it was shown that by introducing a potential for the mimetic non- 
dynamical scalar field, one can mimic nearly any gravitational properties of normal 
matter. It was also shown that in particular, mimetic matter can provide inflaton, 
quintessence and even lead to a bouncing non singular universe. Symmetric static 
vacuum solutions were found for MG in [28]. MG has been extended to f(R) 
gravity in [29, while cosmological viable models derived from such an extension 
are presented in [301 . 

Despite the references above, due to its recent elaboration, MG still presents 
a lack of cosmological and astrophysical application. Our purpose in this work is 
to verify MG at an astrophysical level, by analysing compact objects. By compact 
object in astrophysics we mean a homogeneous distribution of ordinary (or exotic) 
matter with energy-momentum tensor in the form of a perfect non-viscous fluid. 
The key note to start studying the dynamics of such systems is to write the EoM 
for a spherically static metric in Schwarzschild-Droste coordinates, denoted by 
x M = ( t,r,d,ip ). By rewriting the EoM in terms of the hydrodynamic parameters, 
hydrostatic pressure p, energy density p, and mass M (active mass inside the 
star), we will obtain the TOV equations. The matter inside the compact star is 
composed of nuclear matter, sometimes in the form of strange quark matter M 
El|33|, with a density of order p ~ 10 15 pr/cm 3 . It is believed that this matter 
has ordinary statistical properties and can enter a series of second order phase 
transitions m, Eg. 

We will derive TOV equations in the MG framework and solve them for dif¬ 
ferent equations of state (EoS). The paper is organized as follows. In Section [ 2 ] we 
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briefly describe MG. In Section [3] a model of compact star in MG is presented 
from the development of the TOV equations. In Section [4] numerical solutions for 
TOV equations are obtained from both neutron star (NS) and quarks star (QS) 
EoS, while in Section [5] we discuss the results. 


2 Mimetic Gravity 

A generalization of the conformal transformations g^ —¥ fi{x^)g^ was introduced 
by J.D. Bekenstein as disformal transformations [361 : 


9tiv ->■ A' )g llv + B(<j>)Vn<j>V u <j>, X := v (j>. (1) 

Here, <j> is the scalar field and A' the kinetic term. A special case happens when we 
consider B = 0 and A = A( X), which is called MG. It was proven very recently 
that it can be interpreted as the same spirit as the disformal theory ma¬ 
in MG we consider the physical metric g MU as a function of an auxiliary metric 
g M i, and of a scalar field <f>. Therefore one has g^v = (g a ^d a 4>dp<j))g. The MG 
action is therefore 

S = - iq ^ Gn J ^xyJ-gig^^iRlg^ig^, </>)] + Cm], (2) 

for which Cm is the matter lagrangian and Gn denotes the Newtonian gravitational 
constant. Throughout this work we are going to assume units such that c = 1. 

By varying this action with respect to g^ and <j> yields, respectively, the fol¬ 
lowing EoM [251 : 

(G^ - 8nG N T H - (G - 87 rG JV T) 5 /i “/' 3 S a 0^0 = 0, (3) 

-^=d K {V= 9 (G - 8 t iG N T)g KX d x 4>] = 0, (4) 

with G and T being the trace of the Einstein’s tensor and of the energy-momentum 
tensor, respectively. The trace of the Equation © is: 

(1 - g af) da<t>dp<t>){G - 8nG N T) = 0. (5) 

In order to obtain a scenario dynamically different from GR, we suppose that 
(G — 87 t GpfT) 7 ^ 0, and consequently obtain the following equivalent EoM for the 
scalar field: 


g a ^ d a (f>dp(t> = 1. (6) 

This equation is in the Hamilton-Jacobi form and its general solution can be 
obtained by the same method. Generally speaking, cj) can be a complex variable, 
i.e., (j) C. 

It was shown in [26] that such an approach could explain the cold dark matter. 
Here, instead, we will analyze the interior of compact objects from the MG TOV 
equations, which will be derived in the next section. 
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3 Model of a compact star in Mimetic Gravity 

We suppose that the metric of the compact star is static and spherically symmetric, 
with Schwarzschild-Droste coordinates x M = ( t , r, 9, tp) in the following form: 


ds 2 = e 2 ^dt 2 — e 2 X dr 2 — r 2 (dd 2 + sin 2 8 dip 2 ), 


(7) 


for which the coefficients ip and A depend on the coordinate r only. For such a 
metric, the non vanishing components of the Einstein’s tensor read: 



( 8 ) 


0 ) 


G 33 = ~re 2 X [ip' - A' + r(ip" + ip' 2 ) - rip' A'] 
G 44 = — sin 2 8 G 33 , 


( 10 ) 

( 11 ) 


for which primes denote derivation with respect to r. 

By considering the energy-momentum tensor of a perfect fluid, i.e., TK = 
diag(p, —p, —p, — p), with p and p representing the density and pressure of the star, 
respectively, the diagonal components of © are: 



( 12 ) 



(13) 


Also, the trace equation gives us the following: 

2A = 0. 



(14) 


From this differential equation we observe clearly that the auxiliary scalar held <p 
can be a complex function, i.e. cp £ C. We must not be worried about this solution, 
because even in this case, the reparametrization of the physical metric does 
not yield a complex metric. 

Moreover, the hydrostatic (continuity) equation = 0 for v = r yields 


f r =-(p + p)iP'. 


(15) 


This equation does not deviate from the one derived from GR because the auxiliary 
scalar held introduces no extra ghost degree of freedom. The only contribution 
comes from the metric and the huid. 

In order to proceed with TOV equations, we redehne the metric function A in 
terms of a mass function M (r) as follows: 



2 G n M 


( 16 ) 


r 
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From this equation we can calculate the differential change in mass dM, as the 
mass which is stored in a layer with thickness dr in the spacetime, as the following: 


G N dM 

dr 


1 — e _2A (l — 2r\') 


(17) 


with M being the mass of the compact object. We can integrate dM to find the 
total (active) mass inside a region V = {r € [0, /£]}, with R being the radius of the 
object. It was assumed that the mass enclosed inside the compact object is finite 
and a monotonic function of radial coordinate r. 


The plan now is to rewrite Eqs.( 12)-(15 I in terms of { dp/dr,dM/dr,p,p } in 
a dimensionless form. To do so, we make the following sequence of definitions: 
M — > ttlMq, r —> r g r, p -> pM(p/r g , p —> pM@/r g and R —> R/r g , with r g = 
GnMq = 1.47473fcm and Mq is the mass of the Sun. 

Such definitions yield ( |16| )- © to convert to the following: 


^ r dm 
_ m dr 

dr r 2 — 1 
r 


(18) 


Using ( 15|18 ) in (12) yields, for the ( tt ) component: 


2m l~ ^ 1 _ 8t rp 

r 3 ~~ i i- 2 ^’ 


(19) 


while the (rr) and trace equations are, respectively: 


f \ z / f \ ' 1 r dm / 

P \ ( P \ m i- p 

P + Pj \P + p) r 2 2m - 1 p + p 

P 1 _ ,_ 1 _ 8?r(p - 2 p) 

r{p + p) r 3 ^ - 1 2^(1-^) 1 -^ 


( 20 ) 


d<f> _ r g 



We can combine Equations ( |15|19|20|2l] ) and obtain: 


,// _ _ -Ar 3 np - m 

(j)' Am 2 — Amr T r 2 ’ 

i/j" ^ — r 3 + 2 r 2 m) + m(ip'r + 2i/> ,2 r 2 — 1) — ip'r 2 
—ip' 2 r 3 — Aip'r^np — 167r? -3 p = 0, 
m! = — Anpr, 
p = -i>'{p + p), 


( 21 ) 


( 22 ) 

(23) 

(24) 

(25) 


with the following constraints: 2 (f>' 2 m — rr/)' 2 — r g r = 0, r ^ 2m, p + p ^ 0, <j/ ^ 0. 
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An alternative elimination process is to eliminate m first, then p, p, <p and finally 
ip. From this procedure we obtain the following system of differential equations: 


nq 2 / / 3 r> ±f 2 / / 2 2 | if 2 | if 3 r» ±f 2 2 if 
SZpnr (p —Zip r g ip r + tp r g + tp —Zip r g r ip 

—2tp'rgip'r + 2ip'r 2 ip"rg = 0, 

_ 1 0' 3 + (p'rg - 2 rip"rg 

^ 8 nr 2 ip' 3 ’ 

20 0 r ip r g — oip r (f> r g + 50 rip + 70 rip r g 

A if 2 if if 2 I ^ ,/ 3 i// 2 / /2 ^ , /2 3 2 ,/, ff 

—40 r ip 0 + 00 r 0 r g ip — Dtp r r g ip ip 

c // 2 /// Q j/2 2 2 iff oi/2 3 2 iff 

—Z(p rr g (p — 2<p r r g ip — 20 r r g ip 

r\ if2 I /3 3 2 r\if2 2 
— 2<p Ip r r g — 20 Tp 

//4 I j/ 2 3 iff iff r\ ±f2 2 2 / /2 _ 

—zip + Dip r ff r tp tp —Zip r r g tp = 0, 

with the following constraints 


r //3 rr ±f 2 i q j// 2 ( r> // 2 / f2 2 

—50 — 50 r g + 8r0 r g + 20 r g ip r 

+2tp'r^r 2 tp" 4- 2p> r 2 g tp'r — 2tp'r 2 tp" 0, 

^Vo,^o. 


(26) 

(27) 

(28) 
(29) 


(30) 

(31) 


But let us now revert to ( 22]|25 I, which gives us a system of first order autonomous 
differential equations (see |38| for a review of dynamical systems approach in cos¬ 
mology): 


4 > = C, 

/ _ r g — Ar 3 np — m 
k ( 4 m 2 — Amr + r 2 ’ 
ip' = p, 

p' (— r 3 + 2 r 2 m) + m(pr + 2p 2 r 2 — 1) — pr 2 
—p 2 r 3 — Apr^np — 167r?' 3 p = 0, 
m! = — 47rp? , 

P = ~p(p + p). 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


These equations can be written in the form of an autonomous set of first order 
differential equations: 


X' = /(X;r),X = (4>,C,ip,p,m,p). 


(38) 


We use ( 33]|37 ) to find all stationary points of the system. Such points are the 
solutions of /(X c ,ro) = 0, with ro being an arbitrary value of r. For our problem 
we can solve this and find the following stationary point: 


X c — (0 C ,O,0 C , 0,0,0). 


(39) 


We can now linearize the system of equations ( 33||37| ) in the vicinity of this critical 
point. The technique is to make A' = X c + <5X and expand the function / up to the 
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first order of JX. In this case, thanks to the stationary condition, we will obtain 
the following set of first order “linear” differential equations for the perturbations: 

? rt X = JSX. (40) 

dr 

Here, the matrix of the coefficients is called “the Jacobian matrix” and is defined 
by J = -^x\x c - For our case we calculate it explicitly as follows: 


0 1 0 

0 

0 

0 

0 0 0 

0 

r 2 

' 9 

X a r ' 2 

0 

0 0 0 

1 

0 

0 

0 0 0 

—r _1 

-r" 3 

— 16-7T 

0 0 0 

0 

0 

— 4-7T p p r 2 

0 0 0 

—p (o) 

0 

0 


with X c —» 0 and by p p we mean dp/dp. 

To classify the stationary point as stable or unstable, we need to compute all 
eigenvalues of this real valued matrix. An eigenvalue of J is a scalar, such that 
there exists a vector (the corresponding eigenvector) for which JV = AV. The 
characteristic polynomial is given by the following six order algebraic equation: 

A 3 [A 3 r + A 2 — 16 p (0) 7r A r + 4 p (0) 7r p p \ = 0. (42) 


In the above equation we suppose p = p(p) so that p( 0) means p(p = 0). The 
solutions degenerate for A.; = p 2 ,3 = 0, so the system is unstable under infinitesimal 
perturbations X —> X + <5X. Another real valued root is given by the following: 


_ 1 C 4 + 48p(0)?rr 2 + 1 - C 2 
Xr ~3 - C^r - 

with C defined as 


(43) 


C 6 = —72 p (0) 7t r 2 — 54 p (0) n p p r 2 — 1 + 6 Br, (44) 


and 

o2 

— — j-t—= 1024 p(0) 2 n 2 r 4 + 16 p(0) ir r 2 (45) 

— O p (UJ 7T 

—72 p (0) 7t p p r 2 — 27 p (0) 7r p v 2 r 2 — p p . 


4 Solutions for different star models 


We will solve Eqs. ( |15|19|20|2l| numerically with an auxiliary EoS p = p(p) for 
each of the compact star models below. Before starting to solve such equations 
for QS models, it is worth stressing that QSs are, indeed, likely to exist. In [39j , 
for instance, it is presented a form of probing the existence of compact binaries 
containing at least one QS from the observation of gravitational waves. Once QSs 
existence is probed, the ground state of matter will be understood as that of 
strange matter, i.e., a quark fluid with about equal numbers of up, down and 
strange quarks. 










4.1 QS model with linear EoS 


Let us start solving the system ( |15|19|20|2l| ) by adopting the following quark 
matter EoS 1401: 


p = c(p-4B),c = 0.28,0.98 < B < 1.52. (46) 

Solutions for such a system are given, analytically, by the following elementary 
expressions: 


<(>(r) = -2 r g r + Ci, 
m(r) = -r, 

5 1 


p(r) = 


96-7T r 2 ’ 


p(r) = 
ip(r) = 


\r* 


967T 


R 2 ) ’ 


In 


+ 1 L(1 + c)R 2 — cr 2 J ’ R 


< 2.139. 


(47) 

(48) 

(49) 

(50) 

(51) 


It is noted, for the density function p(r) in the vicinity of r = 0, a singularity for 
such a QS model. Furthermore, the metric function i/j(R) = 0. 


4.2 QS model with exponential EoS 

An empirical EoS is proposed in m as the following: 


P_ 

Pc 


,, -b—\ 

a(l — e pc j. 


(52) 


Here, p c = pc = ^/Mq, and a, b are constants. Low density profiles increase 
significantly the pressure and makes the hardenabilit.^\ of the star. The QS with 
linear EoS p = Ap + B can be considered a special case of (521 when p p c . 


—b — 

Table 1 The parameters of EoS for three different stars: — a(l — e pc ). 


Strange star candidates 

a 

| <?a 

b 

a b 

Her X-l 

0.45526 

| 0.0198 

0.0152 

0.00124 | 

SAX J 1808.4-3658 

1.2242 

| 0.0563 

0.0191 

0.00162 | 

4U 1820-30 

1.1255 

| 0.0408 

0.0243 

0.00163 | 


We list in Table 1 above the parameters for three astrophysical objects. In 
the 4U 1820-30 object, the effects of high pressure are much more important than 
the same effects on the other objects. Also we tested the parameters of EoS (521 


x The depth up to which a star is hardened during a matter collapse process. 
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for several initial conditions. We did not observe any significant difference among 


the models. The schematic form of EoS (521 remains the same. In fact, it can be 


considered as a generalized exponential virial EoS proposed in the literature m- 


pk B T 


= exp 


( r-iOO r/ 7> 
Zjm=2 r ^mP 


(53) 


In the EoS above, the coefficients K m and B n can be computed in terms of the 
Virial coefficients. Also, k B is the Boltzmann constant and T is the temperature. 
However, this similarity is just formal. There is a significant difference between 


our EoS (52 I and exponential Virial EoS (531. It is related to the pressure of the 


background p;, ~ p c . In Kenneth’s EoS (53), because it works mainly for gases with 


exponential virial EoS, the pressure of the background pj, = 0. Meanwhile in our 
compact object model, we have a non zero p h value. Another difference is that we 


constructed EoS (52) numerically from the gravitational field equations, whereas 


exponential virial EoS (531 was made from high-density fluid. 

We solve ( |15|19|20|2lf numerically. The results are shown for Her X-1,4U 1820- 
30 and SAX J 1808.4-3658 in the following. 


4-2.1 Her X-l 


For Her X-l, m ~ 0.88, R ~ 7.7 and we expect the pressure to vanish at r ~ R |42l 
[43 [ . Our numerical simulations for dl5|19|20|2lT > are shown in Figs.l and 2. From 
Fig.l, we see that the mass function always increases. Scalar field <fi and metric 
function ip also are monotonic-increasing functions of radial coordinate r. But the 
situation is different for pressure p and energy p, which are monotonic-decreasing. 
At large distances r ~ 8, the pressure p vanishes. 


Mass 



Fig. 1 The graphical plot of mass vs. radius for Her X-l with EoS 



4-2.2 4U 1820-30 


For 4U 1820-30, m ~ 2.25, R ~ 10 and we expect that the pressure vanishes 
for r ~ R M- Our numerical simulations for ( |15|19|20|2l] ) are shown in Figs.3 
and 4. From Fig.3, we observe that the mass function always increases. Scalar 
field <p and metric function ip are also monotonic-increasing functions of the radial 
coordinate r. But the situation is different for pressure p and energy p. Both of 





















10 


IQ - 2 [ P/P 0 






Fig. 2 The graphical normalized plots for pre ssur e p, density profile p, metric function ■?/; and 
scalar field <f> vs. radius for Her X-l with EoS (|52[l. 



these functions are monotonic-decreasing. Moreover, at large distances r ~ 10, the 
pressure p vanishes. 


4.2.3 SAX J 1808-4-3658 


For SAX J 1808.4-3658, m ~ 1.435, R ~ 7.07 and we expect that the pressure 
vanishes for r ~ R 143]. Our numerical simulations for ( |15|19|20|2l) are shown 
in Figs.5 and 6. From Fig.5, we observe that the mass function always increases. 
Scalar held <j> and metric function ip are also monotonic-increasing functions of 
radial coordinate r, as for the other objects which are monotonic-decreasing. At 
large distances r ~ 10, the pressure p vanishes. 


4.3 NS with polytrope EoS 

A polytrope EoS for NS was proposed in the following form m : 

p = Kp 1+ ±, (54) 


with K and F being constants. 






11 



Mass 




Fig. 6 The graphical normalized plots of pressure p , den sity profile p, metric function ip and 
scalar field <j> vs. radius for SAX J 1808.4-3658 with EoS ( ]52[ ). 
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Mass 

M 0 


0.0025 - 



Fig. 7 The graphical plot of mass vs. radius for the NS with poly tropic EoS. 


If we integrate the first law of thermodynamics, we obtain: 


p=p+ „ 


(!) 


(55) 


By setting K = 1 and F = 2, this EoS leads to compact objects with accepted mass 
and radius for NSs 05]• 1 11 fact, it has been recently used in Horndeski’s models of 
NSs [IB]. Our numerical simulations for ( 15|19|20|21 ) are shown in Figs.7 and 8. 
From Fig. 7, we observe that the mass function increases at some points and remain 
constant at others. Therefore mass with polytropic EoS is not a monotonic increas¬ 
ing function. Scalar field <p and metric function ip are not monotonic-increasing 
functions of radial coordinate r. For scalar field, the function is linearly dependent 
on r. The metric function ip is increasing for r < 2.5, monotonic-decreasing for 
2.5 < r < 5.5 and it reaches a local minimum at r ~ 5.5. It is monotonic-increasing 
for r ~ (5.5 — 9.5) and at the radius r ~ 9.5 it becomes discontinuous. In Fig.8 
we draw the graphics for density and pressure. The density p is not a monotonic 
increasing function. It has two minima at r ~ 5,11.75 and vanishes at distances 
larger than r > 17. This would, of course, create a periodic behavior, which would 
be determinable for the pressure and density of any particular point (say the cen¬ 
ter) of the star, but the pressure p would be smaller for a point near the star 
radius and greater for a point farther from the star. In fact, for r ~ 8 the pressure 
p vanishes. 


5 Conclusion 

We derived, in the MG framework, the dynamical equations for a static spherically 
symmetric spacetime, which corresponds to a compact star. A dynamical analysis 
showed that the system is locally unstable near the stationary point. The mass, 
pressure, density, metric function and scalar field were constructed numerically. An 
exact solution for a QS model with linear EoS was discovered. A NS model was 
proposed by assuming a polytropic EoS. The qualitative behavior of the functions 
were studied. Our work exposes new features of MG as an alternative for GR and 
other types of modified gravity theories. 
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p 





scalar field (f> vs. radius for the NS with polytropic EoS. The curves for p and p are periodic, and 
cut the r-axis at the point r ~ 2.5. The maximum value of the pressure is pmax ~ 2.5 X 10 —5 
and of the density is pmax ~ 8 x 10 —4 . 
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